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Abstract 

The overdamped thermal ratchet driven by an external (Orstein-Uhlenbeck) noise 
is revisited. The ratchet we consider is unbounded in space and not necessarily peri- 
odic . We briefly discuss the conditions under which current is obtained by analyzing 
the corresponding Fokker-Planck equation and its lack of stationary states. Next, 
two examples in more than one dimension and related to biological systems are pre- 
sented. First, a two-dimensional model of a "kinesin protein" on a "microtubule" 
is analyzed and, second, we suggest that a ratchet mechanism may be behind the 
folding of proteins; the latter is elaborated with a multidimensional ratchet model. 
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1 Introduction 



The ratchet mechanism to produce directional current or motion has become a 
paradigm of the interplay of nonlinear phenomena with Brownian motion and 
external noise in mesoscopic systems [1,2, 3]. We shall use the term "ratchet" 
for any conservative mechanical system that has an ingrained spatial asymme- 
try; a one-dimensional example is a particle in a saw-tooth like potential. Due 
to the potentiality of these systems to describe mesoscopic systems and their 
processes, such as biological ones [4], one should consider the ratchet to be 
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immersed in a thermal bath, giving rise to Brownian motion. This adds dissi- 
pative and stochastic thermal forces that are related to each other through the 
fluctuation-dissipation theorem[5,6]. As it was clearly pointed out by Feynman 
in his Lectures[7], the Second Law of Thermodynamics implies that a ratchet 
under the above conditions cannot generate directional motion. It is by now 
very well established that in order to obtain current or directional motion, 
it is indispensable to have the presence of an external unbiased time depen- 
dent force. This force can be deterministic, a so-called "rocking" ratchet, see 
Ref. [1], or stochastic in nature, see Ref. [8]. In any case, any of those forces 
should be of zero average in time in order not to produce any biased current; 
in the present article we shall be concerned with external stochastic forces 
only. 

It is the opinion of the authors that the origin of the appearance of the current 
is not fully understood yet, although it is clear what the conditions are for 
its existence. In particular, a lot of attention has been devoted to the case 
in which the ratchet variable is periodic, such as an angle 6 G [0, 2n], and 
many analytic and numerical results have been obtained [2, 9], yielding a clear 
picture of how current occurs in a stationary state. We believe the studies 
are not as exhaustive for the cases where the ratchet variable is not periodic, 
namely when the particle can move in all space, namely, when the position of 
the particle can take all real values, x G (—00,00). In this case, even in the 
absence of external forces, the system does not reach an stationary state[5]. 
One of the purposes of this article is to add to the understanding of the 
appearance of the current in the latter situation, by studying a novel non- 
periodic one-dimensional overdamped ratchet; this is done in Section 2. 

Having established the conditions for the production of current in one di- 
mension, we then proceed to the second purpose of the article and extend 
the model to more than one dimension in order to study two biology-related 
problems. One is the case of the motion of a kinesin protein on a microtubule, 
in Section 3, and the other is the process of protein folding in Section 4. The 
first has been studied already by several authors[10,ll] but here we consider 
a protein in two dimensions and the microtubule as a quasi-one dimensional 
structure. Our model of protein folding is at present speculative, and the idea 
is centered on the hypothesis that the energy landscape of the protein has a 
ratchet-like structure; by further assuming that there is an external agent that 
consumes energy, which could be the presence of chaperone proteins[12], one 
has all the ingredients to obtain directed motion. We shall show that, under- 
standing by protein folding the search and finding of a target point or small 
section in the energy landscape, the ratchet mechanism is extremely fast with 
an almost perfect efficiency An interesting consequence of this process is that 
the landscape need not have a funnel structure[13,14]. 
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2 Current in a one-dimensional overdamped ratchet 



In this section we revisit the one-dimensional overdamped ratchet system, con- 
sidering both periodic and non-periodic potentials. The main purpose here is 
to argue that the appearance of a current is a generic property of bounded 
asymmetric conservative systems, in interaction with a thermal bath, due to 
the presence of external (stochastic or deterministic) non-biased time depen- 
dent forces. That is, by showing that the absence of current without the ex- 
ternal force must follow from the Second Law, we argue that there is nothing 
to prevent a current once the conditions required by the Second Law have 
been relaxed. The details of the numerical techniques we use can be found in 
Refs.[15,16]. 



2.1 No current in the absence of external time- dependent forces 



We initiate by writing down the the Langevin equation in the overdamped case 
for a particle subjected to a conservative potential V(x), without an external 
time dependent force, 

dx dV „ , . . 



We consider a ratchet potential with a "quenched" random disorder both in 
the potential barriers and in the spatial periods, namely, 



V(x) = -V^ 



(0 



. / 2tt(x - y. 
sin - 




if yi < x < y i+1 , (2) 

with i taking all integer values and where the positions are chosen at ran- 
dom; y is arbitrarily set equal to 0. Clearly, the spatial periods are given by 
Aj = yi + i —y%. A typical form of the ratchet potential V(x) is shown in Fig. 1. 
Note that this potential is not periodic. The usual periodic ratchet potential 
is V W = V and \i = \ for all i[2\. 

Returning to the Langevin equation (1), 7 is the friction coefficient and f(t) 
is the thermal force exerted by the thermal bath. The stochastic properties of 
this force are 

(/(*)> = and (f{t)f{t')) = 2 1 kT5{t-t'). (3) 
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Fig. 1. Ratchet potential with quenched random disorder both in the potential 
barriers and in the spatial periods. See Eq.(2). 

It is further assumed that this process is Gaussian, so that Eqs.(3) suffice 
to completely specify the process f(t). In the second equation T is the tem- 
perature of the bath and the expression for the correlation of the force is the 
Fluctuation-Dissipation relation for this problem; it imposes detailed balance [5,6]. 

Since the process x(t) representing the position of the particle is Markovian, 
then, knowledge of the conditional probability distribution p(x, t\xo) is enough 
to determine the process. It can be shown that this function obeys the following 
Fokker-Planck equation [5], 



supplemented by the initial condition p(x, 0\x ) = 5(x — x ). 

In this discussion we are considering that the particle can move in all space, 
namely, — oo < x < oo. Thus, the normalization of the distribution, 




(4) 



oo 




(5) 



— oo 



imposes the boundary conditions 



p(x, t) — >0 for x — > ±oo. 



(6) 
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The Fokker-Planck equation, Eq.(4), is a continuity equation for the probabil- 
ity distribution. Therefore, one can read off the probability density current, 



./x IdV , , kTdp(x,t) 



We may define the "total" current as, 



J(t)= J j(x,t)dx 



=--<-> 

= |wo>, (8) 

where the last two lines follow from the definition of the density current and 
the Fokker-Planck equation, Eqs.(7) and (4). 

If there exists a stationary distribution, p s (x), then detailed balance means 
that the stationary density current vanishes, j s (x) = 0, and therefore that the 
stationary distribution is 

Ps ~ e-™*\ (9) 



with (3 = 1/kT. Further, if p(x,t) is not the stationary distribution, then, 
detailed balance implies that as t — > oo, such a distribution approaches the 
stationary one. This is equivalent to the if-theorem. 

In the context of ratchet-like potentials, one is faced with potentials that are 
bounded, that is, V 1 < V(x) < V 2 , for all x. Thus, strictly speaking, a station- 
ary distribution cannot be reached[5]. However, it is clear that stationarity is 
achieved in the following sense, 

j im ^ 1 *I = -/J(V(x)-V (*')). (1Q) 
t-«x>p(x',t) V ' 



Therefore, it must be true that, in the above sense, the probability density 
current j(x, t) and the total current J(t) must vanish as t — > oo. However, this 
does not suffice to prevent an arbitrary current or displacement of the particle, 
as is routinely asserted in essentially all articles dealing with this problem. The 
main point that we want to emphasize here is that, since detailed balance is 
in accordance with the Second Law of Thermodynamics, the total current 
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J(t) must be bounded for all times. Using Eq.(8), this can be more precisely 
written as the following requirement, 



\{x(t)) - x \ 



< A(ar ), 



;n) 



where Xq is the initial position of the particle and X(xq) is the distance between 
the adjacent local maxima of the potential where the particle was initially at 
x = xq. We have exhaustively numerically verified that this holds for ratchet 
potentials such as those ones in Fig. 1, and the results are exemplified in Fig. 
2. 
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Fig. 2. Average position of the particle (x(t)) as a function of time, for the case of 
no external time-dependent force, see Eq.(l). (a) For a disordered ratchet potential, 
(b) For a periodic ratchet potential. In both cases the initial condition is x = —0.04. 
Note that the particle, on the average, never leaves the well where it started, see 
Fig. 1. 

For ratchet potentials this means that on the average the particle cannot 
leave the local well where it was initially; if it did, nothing would prevent 
the particle to "jump" to another well, and so on, thus moving an arbitrary 
distance. In other words, Eq.(ll) expresses that fact that it is not possible 
to obtain a current that, on the average, could generate motion that would 
transport the particle farther than its initial well. Otherwise, this would yield 
the possibility to perform work on some external load or agent [17] violating 
Kelvin's statement of the Second Law: "A transformation whose only final 
result is to transform into work heat extracted from a source which is at the 
same temperature throughout is impossible" [18] . In the opinion of the authors, 
this is what Feynman, in his Lectures[7], wants to convey in his discussion of 
a ratchet engine. 
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It is important to stress that we consider high enough temperatures so that 
the particle, even though it does not move on the average, it does perform 
normal diffusion. This is shown in Fig. 3. 
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Fig. 3. Second moment of the particle position distribution Ax 2 = (x 2 (t)) — (x(t)) 2 
as a function of time, for the case of no external time-dependent force, see Eq.(l). (a) 
For a disordered ratchet potential, (b) For a periodic ratchet potential. As discussed 
in the text the particle performs normal diffusion. 



2.2 Current in the presence of an external time- dependent force 



In the presence of a time-dependent external force F(t), the Langevin equation 
reads, 

dx dV „ , . . . . . 

^ = -& +/ « +F «- (12) 

where the thermal force f(t) obeys the same properties as before, see Eqs.(3). 
As an example, the external force F(t) may be generated by an Orstein- 
Uhlenbeck process 

dF 

ro^ = -F + CW, (13) 
with ((t) a given Gaussian, white, stochastic process: 

<C(f)> = o <C(OC(0> = /o*(*-0- (14) 



The parameters r and f determine the process. The first one is the correlation 
time for F(t) and fo is a measure of its strength. In order to avoid any bias 
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by this force, one must consider the initial condition F(0) = 0. 



Equations (12), (3), (13) and (14) are equivalent to the following bivariate 
Fokker- Planck equation [9], 



at 7 ox \ ox ox J 

+ 7 w( F+ iM w(x ' F ' t] ' (15) 

where W(x, F, t\xo, 0) is the conditional probability distribution to find values 
x and F for the corresponding stochastic variables, given that at t = 0, x = x 
and F = 0. Initially, it is W(x, F, 0\x , 0) = S(x - x )S(F). 

By simple inspection we find that, even if a stationary distribution W s (x,F) 
exists, equation (15) does not obey detailed balance. The "offending" term is 

l£-(-FW{x,F,t)). (16) 
7 ox 



This, of course, only means that F(t) is "external". Namely, F acts on x but 
not the opposite, Thus, there is no mechanism to establish equilibrium among 
the particle, the thermal bath and the external agent. Therefore, the system 
can withdraw energy from the external source and generate motion. In other 
words, there is nothing to prevent the total current from taking any value 
different from zero: 

1 00 °° P)V( \ 

J(t) = — I dx J dF%F,i)^^0, (17) 

^ — oo — oo 



or 

J(t) = j t (x(t))^0. (18) 

In other words, in this case the Second Law does not impose any restriction 
on the appearance of a current. The Second Law will now impose restrictions 
on the efficiency of the process but that is not discussed here. Of course, if 
the potential V(x) does have a "left-right" symmetry, the current vanishes as 
one should expect. 

In Fig. 4 we exemplify the appearance of a current, both for the disordered 
potential and for the periodic one. On top of the current, the particle also 
performs normal diffusion, i.e. for long times < x 2 (i) > — < x(t) > 2 ~ Dt. 
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Fig. 4. Average position {x(t)) as a function of time in the presence of an external 
time-dependent force, (a) For a disordered ratchet potential, (b) For a periodic 
ratchet potential. 

3 A two-dimensional "kinesin on a microtubule" model 

During the last few years there has been an increasing interest in studying 
the statistical behavior of the transport phenomena inside the cell carried out 
by protein motors and it has been proposed by several authors that ratchet 
models may be relevant in the description of these processes[l,2,8,9,10,19]. 
Protein motors are responsible for carrying diverse kind of vesicles from one 
site to another in the cell in a much more efficient way than the obtained by 
simple diffusion[12]. To accomplish this task they consume chemical energy, 
usually stored in the form of adenosine-triphosphate (ATP), to convert it into 
mechanical motion [2] . Due to their dimensions, the erratic collisions with the 
solvent molecules represent a non negligible contribution to the protein motors 
dynamics and appear in the form of friction and thermal stochastic forces. This 
is why they are also called Brownian motors or molecular motors [4]. 

Molecular motors, as kinesin and dynein, present directed motion only when 
they are attached to microtubules, otherwise they perform standard Brow- 
nian motion. It has been also observed that they do not hydrolyze ATP at 
an appreciable rate unless they are attached[2]. Experimental observations 
of kinesin and dynein motion on microtubules [20] have shown that kinesins 
move mainly along one single protofilament, while dyneins visit often several 
protofilaments. One important characteristic of the microtubules is that they 
have an intrinsically periodic and asymmetric structure[12] . 

Most of the models motivated by these experimental results were initially 
restricted to the one dimensional case but it has become clear that more di- 
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mensions are needed[ll,21,22,23,24,25,26,27,28,29]. Some of these works [28, 29] 
approach the problem imposing ad hoc asymmetric probabilities thus obtain- 
ing transport. The others use Langevin equations and, in particular, those 
in Refs. [11,21,22,23] have introduced more detailed models to understand the 
particular behavior of kinesin in microtubules. 

As an step forward in the description of these processes, we are considering 
here a model for a kinesin on a microtubule as a ratchet in a two dimensional 
space, immersed in a thermal dissipative bath, and subjected to an external 
force of zero mean varying stochastically in time. At present we study general 
statistical properties of the model; a more detailed study and its comparison 
with real kinesin proteins will be presented elsewhere. Nevertheless, the results 
shown here resemble qualitatively what is observed in actual experiments with 
motor proteins and microtubules [12]. 

In order to model the interaction forces between the motor protein and the 
microtubule a two-dimensional potential is proposed, whose shape resembles 
a long attractive filament (or channel) with a typical ratchet structure inside 
on the longitudinal direction, as shown in Fig. 5. The potential V(x, y), called 




Fig. 5. Two-dimensional potential with an infinite attractive filament (or channel) 
along the x coordinate and a ratchet structure inside of it. The attractive filament 
has a mean depth of V\ = 15 and a width determined by a = 2, while A = 1 
represents the period of the ratchet and Vq = 2 its amplitude. The longitudinal 
projection of this potential is similar as the used in one dimension models. In the 
transversal direction, — oo < y < oo. 

for simplicity the filament or microtubule potential, is given by the following 
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function 



V(x,y) = -e 



-y 2 /- 2 



+ V (sm (— ) + - sm (— ) + - sm (— )) j , (19) 



where V\ is the mean depth of the attractive filament and its width is deter- 
mined by a; A represents the period of the ratchet and V® the amplitude. The 
system is immersed in a thermal bath at temperature T and there is present 
an external random force with zero mean that represents in the model the 
consumption of ATP by the protein motor. Depending on the temperature of 
the bath and on the statistical properties of the external force, the particle 
may be inside or outside the microtubule, i.e. \y\ < a or \y\ > a . Inside, it 
will feel the effect of the ratchet potential and may generate directed motion; 
outside, it will perform free Brownian motion. 

The dynamics of the model are represented by the following coupled Langevin 
equations for a Brownian particle moving in an x — y coordinate system: 

•£ = -^ + (20) 

and 

7 |_av^) +/>(i) + m (21) 



where 7 is the friction coefficient and fi(t), with i — x,y, is the force exerted 
by the thermal bath. The stochastic properties of this force are 

(/i(t))=0 and (f i (t)f j (t'))=2 1 kT8(t-t')5 ij . (22) 



It is further assumed that these processes are Gaussian. The external force is 
Fi(t) with % — x,y. Due to the lack of information with respect to the actual 
statistical properties of ATP consumption by the kinesin, and which is one 
the main difficulties to make precise predictions, we have decided to employ a 
stochastic force generated by an Ornstein-Uhlenbeck process 

r ^- = -F i (t) + Ci(f), (23) 



with Q(t) a given Gaussian, white, stochastic process 

(&(*)> = and (Q( t )Q(t')} = f 5(t-t')5 ir (24) 
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The parameters tq and fo determine the process and represent the correlation 
time and the magnitude of the external force, respectively. As it can be ob- 
served, the relations for fi(t) and Fi(t) are exactly the same for both directions 
and they are ruled by the same parameters. In consequence, they act isotrop- 
ically. However, each of them arise from an independent Gaussian noise, and 
thus, there is no correlation among them. 

We now present the main statistical properties of this model. In all the re- 
sults presented here the initial conditions are x(0) = y(0) = 0. Variations of 
the different parameters, a thorough description of the behavior of the first 
and second moments of the position distribution of the particle, as well as a 
comparison with actual motor proteins, cannot be presented here due to the 
briefness of this report and will be discussed elsewhere. 

As one should expect, the presence of the external time dependent force pro- 
duces a net current along the x direction due to the asymmetric properties 
of the ratchet inside the filament. In the present case, the total current is a 
two-dimensional vector, 

J(t) = J x (t)i + J y (t)j = ±(x(t))i + l(y(t))j. (25) 



In Fig. 6 we show the average x-position of the particle as a function of time, 
the current J x (t) being the derivative of such a curve. The current along y is 
zero, J y (t) = 0, as expected. The current occurs, however, in an interesting 
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Fig. 6. Average x-position of the particle as a function of time. 

manner. First, we realize that the particle feels the effect of the ratchet only 
when it is inside the filament. Outside, it performs diffusive standard Brownian 
motion. The net result is that the current tends, asymptotically and very 
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slowly, towards zero. However, for any finite amount of time, the current is 
different from zero. Thus, a net transport is always realized. Another form 
of saying this is that, even though the particle may stay a long time outside 
the filament, it effectively has directional motion since it eventually returns 
to it. This is better seen in Fig. 7, where we present four "snapshots" of the 
position probability distribution p(x, y, t) for different times, and obtained 
from 2000 realizations of the process. The distribution acquires an arrow-like 
structure because of the ratchet within the filament and the centroid moves 
always to larger values of x, i.e. the distribution appears to be dragged by 
the filament. The shape of the position distribution p(x,y,t) suggests also a 
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Fig. 7. Snapshots of the position probability distribution p(x,y,t) for four different 
times, obtained from 2000 realizations of the process. 

peculiar behavior of its second moment, 



Ar 2 = Ax 2 + Ay 2 

= ((x 2 ) - (x) 2 ) + ((y 2 ) - (y) 2 ) . (26) 

In Fig. 8 we show these three deviations. The total second moment Ar 2 follows 
a power law t a with a ~ 1.5 after t ~ 100 time units. That is, the system effec- 
tively presents superdiffusion, in contrast to the one-dimensional case where 
it is always normal diffusion. The super diffusive behavior is also present in 
the rc-direction, while in the y direction the diffusion is normal. 
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Fig. 8. Total second moment Ar 2 and its components in x and y, see Eq.(25), as a 
function of time. Note that both Ar 2 and Ax 2 show super diffusive behavior, while 
Ay 2 shows normal diffusion. 

The explanation for this behavior appears to be simple: When the particle 
is bound to the filament it feels the effect of the ratchet and tends to move 
preferentially in the positive x direction. When it is outside of the channel 
performs a normal random walk, with normal diffusion, but it does not move 
preferentially to any direction. Thus, for those times when it is outside the 
filament, even though appears statistically to be dragged by the filament as 
mentioned above, the particle lags behind the centroid, hence effectively pro- 
ducing a wider distribution and a larger second moment. It is interesting to 
point out that the superdiffusive behavior is a consequence of "standard" 
noises coupled to a highly nonlinear process and does not arise from an exotic 
stochastic process[30]. 



4 A ratchet mechanism behind protein folding? 

Following the classical work of Anfinsen[31] on protein folding, Levinthal[32] 
argued that if the multidimensional (free) energy landscape had a "golf course" 
like shape, with the hole being the protein native state, it would take essen- 
tially forever for a protein to fold correctly, provided the search were per- 
formed at random, namely, by Brownian motion. Since proteins in vivo fold 
extremely fast and emciently[12] the concept of a funneled energy landscape 
was developed [13, 14], with the further possibility of folding pathways in which 
the protein folds following a path towards or within the funnel leading to the 
native state. But very importantly, the process being driven essentially by free 
energy differences, that is, by "falling" into the state of lowest free energy. This 
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idea has been central in the study of protein folding and we do not pretend 
to review all the advances made by the many groups working on this field [33] 
nor to point out any possible pitfalls of the theories. Rather, we would like to 
point out an alternative way out to the so-called Levinthal paradox, without 
necessarily requiring a funnel leading to the native state. We shall appeal to a 
ratchet mechanism for the process of protein folding. This model requires, in 
addition to having a ratchet-like energy landscape as we discuss below, that 
the folding is a non-equilibrium process being driven by an "external" source. 




Fig. 9. Two-dimensional "free energy" landscape of a protein, see Eq.(27). The 
landscape is flat on the average with a ratchet-like structure on top of it. The deep 
well in the center is the "native" structure of the protein. 

The model consists of a "particle" moving in a multidimensional space, the 
coordinates representing the different configurations of a protein, or effective 
degrees of freedom[36]. The particle is affected by a conservative ratchet poten- 
tial which maybe flat on the average (i.e. no need of a funnel) and with a "small 
and deep enough hole" representing the native state. On top of the average po- 
tential there is an asymmetry that makes the potential toward the native state, 
different than in the opposite direction. Fig. 9 is a two-dimensional example 
for visualization, but we have realized many multidimensional potentials with 
the same type of "toward-away" asymmetry, see below. Although we cannot 
justify the existence of such a potential, one may argue that the regularity 
specially in the secondary structures, at least at first order, suggest a kind of 
"universal" regular potential. Our proposal, undoubtedly not quite justified, 
is that such a regularity may have a ratchet-like form; after all, any coiling 
structure has an asymmetry since it distinguishes between "right" and "left" . 
The presence of the thermal bath needs no justification and simply represents 
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the effect of a viscous environment at a fixed temperature. With only the last 
two ingredients, the ratchet potential and the bath, the "particle" will obey 
Levinthal paradox and would not find its way in a reasonable time, specially in 
a multidimensional space. However, if the there is an external source of energy, 
with random properties and zero bias, acting on the particle, the latter will 
definitely find its way toward the native state and, on the average, could do so 
in a short time. The efficiency of the process, of course, will depend on the val- 
ues of the different parameters. What is the origin of the external source? We 
may generically argue that processes in living organisms are through states of 
non-equilibrium with all sorts of gradients of different physical properties, and 
that the maintenance of those gradients may be traced back to the consump- 
tion of ATP. In a more specific way, although it is not clear that occurs in all 
proteins, it is known that many proteins fold assisted by chaperone proteins, 
which in turn, consume ATP [34]. Thus, the proposed mechanism requires the 
presence of an external source, which in our opinion, is physically appealing 
since Life needs such a source for occurring as required by the Second Law. 
Nonetheless, we cannot say exactly what the mechanism is at the level of the 
protein. 

The mathematical model is, therefore, a multidimensional version of Eqs.(20) 
and (21) for the kinesin-on-a-microtubule model, but now the ratchet potential 
has an overall ( topologically equivalent) spherical symmetry, with a "towards- 
away" asymmetry. For the numerical results we present below we consider 
periodic potentials such as the following, 

V(xi, X2, ■ ■ ■ , Xd) = Vo [sin(27rr/A) + a sin(47rr/A) + 6sin(67rr/A)] 

-V x e' CT \ (27) 

where d is the dimension of the space, namely, the number of effective degrees 

of freedom of the protein; r = (x\ + x\ H V x 2 ,) 1 ^ 2 ; and a, b , c, V and V\ 

real positive numbers. The last term in Eq.(27) represents the "hole" of the 
native state. We want to insist that the ratchet potential need not be peri- 
odic, that is, the saw-tooth structure can have "quenched" random distances 
between maxima and minima, such as in the one-dimensional case, see Fig. 
2. Furthermore, the landscape does not have to be flat, namely, it can have 
valleys and peaks, wider and higher or lower than the average peaks of the 
ratchet, and the particle can still climb over the peaks provided they are not 
above the stall load of the ratchet [17]. Full analysis of these cases are beyond 
the briefness of this report and it will be presented elsewhere. 

Fig. 10 shows a typical example of the many different cases we have analyzed. 
It corresponds to a "protein" with 6 degrees of freedom. We show the average 
evolution of the protein coordinates Xi(t), with i — 1, 2, . . . , 6 over 100 runs for 
the same initial condition. In the presence of the external force, and for those 
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parameters, we found more than 95% "folding" , namely, the particle reached 
the "native" state in the arbitrary maximum time of 15000 units. The average 
time of folding can be read off the graph and corresponds to 4000 unit times 
approximately, very short even in CPU time scales. When the external field 
was turned off, needles to say, we found 0% "folding", and we believe it will 
essentially never find its way no matter how long we run the program. 




Fig. 10. Average particle position (xi(t)) in a six-dimensional space, i = 1,...,6, 
as a function of time for 100 realizations and same initial conditions. The spatial 
period of the ratchet is A = 1. (a) In the presence of the external time dependent 
force, the "protein" folds, (b) In the absence of the external time dependent force, 
the "protein" essentially never folds. 
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5 Remarks 



In this article we have revisited overdamped ratchet systems in one dimension 
and extended it to several dimensions in order to study biology related sys- 
tems. For the case of one dimension we analyzed a novel disordered ratchet for 
the purpose of making explicit the fact that the space in which the systems 
move do not correspond to a periodic variable, as it is routinely done in the 
literature. For these systems, strictly speaking, do not exist a stationary state. 
The main result that we want to stress here is the fact that, in the absence 
of external time dependent forces, stochastic or deterministic, there cannot 
be a current for all times and not only as an asymptotic condition. This is 
a requirement of the Second Law. Once there is an external time dependent 
force, and the system shows an intrinsic asymmetry, such as a "left-right" one, 
there is nothing to prevent a current. The Second Law now acts by setting 
limits for the amount of work released, namely on the efficiency of the process. 
We did not address the latter limits although several workers have discussed 
this point [17]. 

The ratchet mechanism has become a potential candidate for many biological 
process, since, on the one hand resembles a thermal engine at a mesoscopic 
level, but on the other it does not appear to need a specific design such as 
the man-made engines. That is, by its inherent capacity to deliver work in 
the form of directional motion, one may be prompted to use it as a model in 
somewhat obvious biological situations in which a transport process, of some 
kind, is present. This has been the case of the motion of kinesin proteins on 
microtubules, where the protein appears to truly "walk" using ATP as the 
necessary input energy; there are other cases where ratchets have been used 
as models of intracellular motion, or as molecular motors, such as in the mito- 
sis of the cell [35]. Granted, all the models are still at a very premature level of 
comparison with actual situations, mainly because the complicated biochem- 
ical processes involved in real life can hardly be thought to be described by 
systems with one or two degrees of freedom. Nevertheless, we believe it is 
worthwhile to keep exploring these models not only by its potential relevance 
in biological systems but for their own sake. Here, we have presented a two 
dimensional model of a kinesin in the presence of a quasi one-dimensional 
microtubule. The typical walks of the particle, wandering around until they 
hit the microtubule and then directionally moving, show a striking similarity 
with actual kinesin motion, as seen in videos of these systems[37]. Although 
one can always try to fit these simple results to known conditions, such as 
speed of kinesin walks or the energy released by the consumption of ATP, 
we believe further research is needed in two aspects of the model: first, we 
must confidently know how to relate the ratchet potential with the actual pe- 
riodic structure of the microtubule and the form in which the protein attaches 
to it, and second, how to describe the consumption of ATP by means of an 



18 



appropriate statistical process. 

The third section of this article on considering protein folding as being driven 
by a ratchet mechanism is highly speculative, but we believe it is worthwhile to 
pursue it because, again, protein folding appears as a process with directional 
motion. Traditionally, this process has been thought to be driven by free energy 
differences leading towards a minimum, similarly to a chemical reaction[13], 
and thus the idea of a funnel in the free energy landscape of the protein. The 
presence of catalyzers and/or chaperone proteins may accelerate this process 
and their presence is welcome in the theory. Here, by following the idea that 
Life processes consume energy to yield their products, similarly to a man- 
made thermal engine, we have speculated that protein folding is driven by an 
"engine" in which the protein itself is part of it. This now makes a necessity 
the presence of choperones and/or catalyzers that in turn consume energy. At 
this moment, the hardest part to justify in our model is the ratchet structure 
of the energy landscape, and the best that we can say is that it is motivating 
to see that the secondary structure of proteins is universally made of a-helices 
and /3-sheets, structures with certain periodicity and asymmetry. On the other 
hand, we do not have evidence neither pro nor con that a ratchet structure is 
present since this would have to be seen in the highly multidimensional energy 
landscape. As mentioned before, the "folding" with this mechanism does not 
need a funnel, but of course, the interplay of both properties would make the 
process even more efficient. The simulations that we have performed, so far up 
to a landscape in 6 dimensions, are extremely encouraging since they can easily 
be tuned to a very fast and almost 100 % "folding". Extending to an arbitrary 
number of dimensions only requires longer computer time and, at present, 
does not add anything fundamentally different. What is really needed is the 
actual form of an energy landscape, a hard problem with many researchers 
very much interested in it. However, without yet knowing the actual form of 
an energy landscape in its many dimensions, we can make models in which the 
ratchet is not periodic and in which the landscape have valleys and peaks to 
test the efficiency of the search for the "native" state. We shall present those 
results in future contributions. 
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